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Abstract. We propose to use a multimode optical fiber with a D-shaped cross 
section as a privileged system to image wavefunctions of a chaotic system. Scar 
modes are in particular the subject of our investigations. We study their imprints 
on the statistics of intensity and we show how the introduction of a localized gain 
region in the fiber is used to perform a selective excitation of scar modes. 

1 Introduction 

Experimental wave chaos has naturally emerged from the demand to achieve simple experiments 
to study quantum chaos. Since 1980, a lot of experiments have been performed in different sys- 
tems and with different natures of waves. Some years ago, as we were interested in the study 
of the spatial properties of waves that propagate in a chaotic medium, we asked the follow- 
ing question: how to image wave chaos experimentally? In the litterature, some nice ways to 
vizualize the wavefunctions of chaotic billiards have already been proposed. An exhaustive pre- 
sentation of wave chaos experiments is reported in the Stockmann's book[lJ, we briefly cite the 
main experiments related to wave chaos imaging. Spatial distributions of intensity have first 
been observed in microwave cavities by Sridhar et al [2] showing a selection of wavefunctions 
associated with some particular periodic orbits. By measuring the field displacement of vibrat- 
ing plates, Schaadt et al [3] have investigated the statistical properties of the wavefunctions of 
a Sinai stadium plate. In the afore mentioned studies, the measurements were nondirect: the 
mapping of the wavefunctions was obtained by performing a point-to-point scan of the surface 
(of the microwave cavity, or of the aluminium plates). 

A direct visualisation of the wave patterns of a stadium has been achieved on the surface of 
water enclosed in a stadium-shaped container by Blumel et al\A\. This type of experiment 
prevents any quantitative analysis of the statistical properties of the field, but has produced 
interesting results on the direct observation of scarred patterns [5]. 

Optical waves also constitute an evident candidate to image wave chaos. One of the first exper- 
iment with optical waves invoking the properties of classical chaos to explain the behavior of 
waves has been proposed by the group of A. D. Stone[6]. It has been shown that the lasing prop- 
erties of deformed droplets can be understood via a transition to chaotic ray dynamics. Until 
2002, no direct observation of the wavefunctions of a chaotic optical billiard had been realized. 
With the aim of studying the spatial properties of waves in a chaotic system, we have then pro- 
posed an experiment of wave chaos in an optical fiber. The experiment, simple in its concept, 
offers a direct access to the spatial distribution of intensity by imaging the near-field intensity 
at the fiber output. Unlike previous experiments, this experiment also gives the opportunity 
to investigate the properties of a field being the result of the propagation of a superposition 
of several modes. In this paper, we present an overview of our recent results on the spatial 
properties of the wavefunctions of a chaotic multimode optical fiber. From the experimental 
validation of Berry's conjecture on the generic ergodic behavior of wavefunctions of a chaotic 
billiard [7J to the direct observation of a scar mode, optical fibers provide a good experimental 
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system to image wave chaos. 

This paper is organized as follows. In section [2j the experimental setup is described and its 
relevance in the field of wave chaos is presented. Our main results about the experimental ob- 
servation of the typical behavior of wavefunctions in a chaotic billiard are reported in section 
12.21 In order to investigate the scar wavefunctions, we have developped two approaches: the 
first one is dedicated to the signature of a scar wavefunction on the statistics of the field and 
intensity, the other one to a selective excitation of scar by the way of optical amplification. 
Theses recent studies are described in section [3l 



2 Presentation of the experimental system 
2.1 Description of the system 

The multimode optical fiber designed and fabricated in our lab for our studies is briefly described 
in this section. The transverse cross section of the fiber is a truncated disk fabricated from a 
1-cm-diameter preform of silica that has been cut at half its radius and pulled at a temperature 
low enough to maintain the D-shaped cross section. The final dimensions of the fiber after 
pulling are 120 fim for the diameter and 90 fim for the truncated diameter (see Fig. [I]). These 
dimensions differ from the typical 10-/im- diameter of a single mode optical fiber used in the 
field of optical telecommunication. The index of refraction of the silica core n co at our working 
wavelength is around 1.458. It is surrounded by a black absorbing silicon cladding with an index 
of refraction equal to n c i a d=lA53. The cladding is tinted in order to prevent light from being 
guided during the propagation. 




silicone cladding 

(a) (b) 
Fig. 1. Microscope photography of the D-shaped fiber(a) and characteristical parameters of the core(b). 



As the characteristic dimensions of the fiber are large compared to the optical wavelength, a 
collimated light beam that propagates along the fiber can be approximated, in the geometrical 
limit, to a ray that is reflected at the core/cladding interface along the propagation in the fiber. 
A formal derivation of the eikonal equation from the 3D Helmholtz equation is given in [8] and 
justifies the analogy of our system with a chaotic billiard. Indeed, the evolution of a ray along 
the fiber (in the direction of z) can be associated to a time evolution on the transverse cross 
section of the fiber (see Fig. [2j). 

The so-called D-shaped billiard has been proved to be a fully chaotic billiard by Bunimovich[9] 
and ReefTU]. as soon as the truncated diameter D t is larger than half the radius. If the optical 
fiber can be seen as a wave billiard, it also constitutes an interesting analogous quantum system. 
Indeed, the propagation of the field along the fiber is described by the Helmholtz equation. In 
the weak guidance regime, i.e. for a small difference of indices between the core and the cladding 
[(nco— n dad) / n c i a d «l],a given polarization is nearly perfectly preserved along the propagation 
thus enabling to use the scalar form of the Helmholtz equation: 

(A + d zz )1>(r, z) + n(r) 2 ^(r, *) = (1) 
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Fig. 2. Light propagation along the fiber in the geometrical limit of rays. 



with ko the vacuum wavenumber and n(r) the index of refraction in the medium. In the paraxial 
approximation, this equation reduces to the following form: 

iPMir, z) = [~A + V(r)]0(r, z) (2) 

where (3 co = n co ko and (f>(r,z) is obtained from ^(r, z) = 0(r, z)e % ^ coZ . The paraxial approxi- 
mation is associated to the propagation of light close to the optical axis and is fulfilled as soon 
as the angle of the guided ray is much smaller than the cutoff angle max which is given by 
sin max = yjl - (^dad/rico) 2 - 

Equation [2] is then formally equivalent to a pseudo-time independent Schrodinger equation. As 
a consequence, the D-shaped fiber appears to be an ideal wave chaos paradigm of a quantum 
chaos experiment. 



2.2 Typical results 

One advantage of the experiment with the optical fiber is that the experimental setup is con- 
ceptually simple. A beam from a Helium- Neon laser source at XHe/Ne — 628 nm is collimated 
and enlarged so that its size becomes large enough compared to the transverse size of the fiber 
to be in a quasi-plane wave illumination. The main parameter of the experiment is then the 
angle between the beam with the optical axis at the fiber entrance. Indeed, this angle gov- 
erns the value of the transverse wave number in the fiber K, t by the way of the simple relation 
Kt = &on co sin# and, as a consequence, the modes that will be excited and guided in the fiber. 
At the output of the fiber, one can image either the near-field (NF) intensity using a microscope 
objective or the far-field (FF) intensity in the focal plane of a suitable lens. Typical observa- 
tions associated to a value of the transverse wave number n t — 36/i? are presented in figure 
[U As diffraction occurs due to the finite aperture of the fiber, several transverse wavenumbers 
around the main transverse wavenumber n t contribute so that more than one hundred modes 
are superposed at the output of the fiber. The NF figure [3](a) shows a speckle pattern: the 
intensity is statistically uniformly distributed over the whole cross section of the fiber. This 
figure is associated to the superposition of several modes each with a speckle-like behavior as 
predicted by M. V. Berry [7]. A nice experimental evidence of Berry's conjecture is given by the 
FF intensity figure [31(b). In this figure, the intensity tends to fill a ring, whose mean radius is 
centered on the input transverse wavenumber n t and its width is determined by the number of 
contributing modes. This illustrates that the field is constituted by the random superposition 
of a great number of plane waves each with a fixed transverse wave number but with random 
directions: the field of this generic so-called speckle mode can then be seen as a real Gaussian 
random variable. 

The field that propagates along the fiber can be decomposed on the eigenmodes' basis {0 n }' 

^(r,z) = ^c n n (r)e^ (3) 

n 

where c n is the coupling coefficient between the input field ipo( r ) and the eigenmode n , and j3 n 
is the propagation constant of the mode (j) n ((3 being defined as f3 2 = (3 2 Q — k 2 ). In the paraxial 
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Fig. 3. Near-field intensity (a) and far-field intensity (b) for n t = 36 /R. 



approximation, n t <C /? co , and the exact expression of the decomposition of the field on the 
eigenmode's basis is: ^(r, z) = J2 n c n</>n(r) ex.p(i({3 co - k\I (2(3 co ))z). 

We should mention that our experimental system is an open cavity due to an absorbing cladding. 
It is now well established that open systems support non-real eigenfunctions [TTTT2] , whose 
complexness can be measured by the parameter g 2 , 

9 <(He0 2 » 1 ' 

(with the condition that (9m05Re0) = 0). In Ref. |11|12] . it has been shown that q can be related 
to the loss rate due to the leakage at the boundary. In the case of our fiber, measurements of 
attenuation along the fiber lead to q- values less than 10 -2 . 

Thus, legitimely, the field of each mode can be seen as a real random Gaussian variable, one can 
easily deduce that the distribution of intensity of one mode I = \(j)\ 2 follows a Porter-Thomas 
distribution. By superposing several modes, a Poisson distribution for the intensity of the field 
ip(r,z) is then expected [T3]. The figure [H represents the probability distribution of intensity 
for the superposition of modes of the figure [3] and shows a good agreement with the Poisson 
prediction. 




"0 12345678 
1/ <I> 

Fig. 4. Probabiliy distribution of the intensity associated to the figure [3]J a). The Poisson's distribution 
prediction is plotted in dotted line. 



We have experimentally investigated the series of scar modes associated to the shortest unstable 
periodic orbit that is the 2-bounce periodic orbit along the axis of symetry of the fiber's cross 
section. One can expect to observe a constructive interference effect building upon this periodic 
orbit if the following quantization condition for the value of the transverse wavenumber is 
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fulfilled: 

KtC = 2np + | + 2Aip (5) 

where C is the length of the periodic orbit, p an integer, tt/2 is a phaseshift due to the existence 
of a self focal point on the 2-bounce periodic orbit, and Aip is the phaseshift due to reflection 
at the core/cladding interface. This phaseshift also depends on the value of the angle of the 
beam on the interface so that one has to solve the following transcendental equation to obtain 
the quantized value of the transverse wavenumber: 

K t C = 2tv P + | + 8tt arctan \\ '" co l2 " cLad ( ^ ) -1 (6) 

In the optical context, scar modes along the 2-bounce periodic orbit (PO) can be described as 
the modes of an unstable dielectric Fabry-Perot cavity. 
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(a) 



(b) 



Fig. 5. Near-field intensity (a) and Far- field intensity (b) for a scar mode of order p = 5 with n t 
10.35/JS. 



By solving the Eq. El one obtains quantized values of the transverse wavenumber associated 
to scar eigenmodes along the 2-bounce PO. For very precise illumination conditions, we have 
successfully excited a few modes around the scar mode of order p = 5. The NF and FF intensity 
observed at the output of the fiber are presented in figure [5l 

Even if several modes are excited, the intensity tends to localize along the direction of the 2- 
bounce PO emphasizing the main contribution of the scar mode in the superposition. In the FF 
figure, two symmetric peaks are observed and show that mainly one transverse wave number 
propagates along a well-defined direction, i.e., the direction of the 2-bounce PO. These obser- 
vations have constituted the first experimental evidence of scar modes in an optical fiber [T4j. 



3 Looking for optical scars 

The previous results have revealed the optical fiber as an ideal system to study the spatial be- 
havior of waves in a chaotic billiard. If the properties of the field associated to a superposition of 
modes have been studied and analysed [8], the individual eigenmodes cannot be simply observed 
due to diffraction at the fiber input. In particular, looking for an individual scar eigenmode ap- 
pears to be a hard task. In order to investigate these rare eigenmodes, we have developped two 
approaches. The first one is dedicated to the signature of scar modes on the statistics of the 
field and intensity. Kaplan [15] and later Ishio[l6] have studied the deviation from the Gaussian 
statistics due to the presence of scars. In section [3TT| we present our preliminary results on the 
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imprint of scar on the statistics. The second approach refers to the ability to perform a selec- 
tive excitation of scar by the way of optical amplification. The problem of the introduction of a 
gain region in a system that has special properties of propagation also encounters the concerns 
of recent works on microlaser cavities or on disordered systems that will be made explicit in 
section 13.21 
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Fig. 6. Intensity of the field associated to calculated scar mode (a) and speckle mode (b) and their 
respective probability distributions of intensity (c) and (d). 



Figure [6] displays the intensity of two calculated eigenmodes: the strong localisation of intensity 
along the 2-bounce PO observed for FiglBJa) corresponds to a scar mode, whereas Fig|6jb) is 
associated to a generic speckle mode. The probability distributions of the intensity associated 
to these two characteristic modes compared to the Porter-Thomas distribution deduced from 
the Gaussian variable theory are shown in figure [BJc) et (d). If the probability distribution of 
the intensity for the speckle mode is in good agreement with the Porter-Thomas distribution, 
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a strong deviation is observed for high values of intensity in the case of the scar mode. This 
deviation is related to the 'extra-density' of the intensity along the direction of the periodic 
orbit: high values for scar modes are more probable than for speckle modes. 
As scar modes seem to mark the intensity distribution, one can expect to look for the presence 
of scar modes in a superposition of several modes by analyzing the distribution of the intensity. 
This assumption has been studied in the context of open billiards by H. Ishio et. at. who have 
revealed the influence of 'regular states' on the probability distribution of the superposition. 
Experimentally, a favorable illumination is launched in the fiber in order to excite several modes 
around the scar mode. As the transverse wave number of the fiber ennals 18.6//?. the number 




Fig. 7. NF intensities and their respective probability distributions for an input illumination along the 
2-bounce PO (a), (c) and for the rotated fiber (b), (d). 

Using the calculated eigenmodes, we have artificially constructed a superposition of five modes 
around a scar modes (see Eq|3|) where the coupling coefficients are fixed by the initial field 
^o(r). In one case, the initial field is a plane wave with a transverse wave number associated 
to a scar launched in the direction of the 2-bounce PO. In the other case, the same plane wave 
is launched at an angle of 45 degrees. This launching condition is numerically equivalent to a 
rotation of the fiber around its axis in the experiment. In the latter case, the coupling coeffi- 
cient of the scar mode is very small. In both cases, we consider the field resulting from a 1-m 
propagation. 
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Fig. 8. Intensity of the field associated to a scar (b) and its four neighouring speckle modes (a), (c), 
(d), (e). Probability distribution of intensity for a superposition of these modes (shown in the inset) 
associated to an input illumination in the direction of the scar (f) or an arbitrary input illumination 

(g). 



The figure [8] presents the NF intensity and the probability distribution of the intensity for the 
two different launching conditions. The NF intensity is statistically uniformly distributed over 
the whole transverse cross section showing no difference between the two cases. Nevertheless, 
a deviation for high value of intensity is also observed for the "scar" illumination whereas the 
Poisson distribution is observed for the rotated initial condition. 

As expected, as no scar is excited, no evidence of scar is observed. These results are prelimi- 
nary and will be completed by a more systematic experimental investigation in a forthcoming 
publication. 



3.2 Optical Amplification of Scars 

3.2.1 General context 

Kudrolli et al have studied the influence of nonlinear Faraday waves on the selection of scar 
modes in a stadium shape container [23], but only very few investigations have been done on 
the behavior of the chaotic wavefunctions in a nonlinear system. Recently, chaotic microlasers 
associated to wave chaos in the presence of gain have been the support of a great number of 
investigations [25 26J. These studies have shown that a lasing effect occurs for scar modes thanks 
to their high quality factor, and the privileged directions of emission associated to scar modes 
appear to be interesting to provide highly directional lasers. 

In the context of wave propagation in complex media, some of the spatial properties of waves in a 
disordered system are somewhat similar to those met in chaotic systems. For instance, numerical 
investigations have recently shown that the amplified modes in a random laser are the modes 
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of the passive disordered system [28]. The same authors have shown that the identification 
between lasing modes and passive modes holds even when the gain region is smaller than the 
mode size [29 . 

We have shown that a passive optical fiber with a D-shape constitutes a powerful exper- 
imental system for studying wave chaos. But above all, multimode D-shaped fibers are also 
widely used to achieve optical power amplifiers [31] in which a telecom signal (called the sig- 
nal) is amplified at the cost of an ancillary signal (called the pump). Usually the optical index 
of the Erbium or Ytterbium doped-core is higher than the index of the surrounding D-shape 
cladding. The signal propagates in a single mode doped core, and, by guiding the pump, the 
multimode fiber only plays the role of an energy reservoir. By some appropriate modifications 
of the fabrication process, it is possible to obtain a negligible index mismatch between the small 
active region and the multimode fiber Q, at least sufficiently low to ensure that no guided mode 
exists inside the doped area. Then, the signal and the pump propagate on the modes of the 
entire section of the multimode D-shape fiber. In the following, we present realistic numerical 
simulations of such a chaotic fiber with an Ytterbium-doped region. The simulation is based on 
the Beam Propagation Method we already used successfully for simulating optical amplification 
[32] . We show how localized gain may selectively amplify a family of scar modes. 



3.2.2 Results 

In our multimode fiber, selective excitation of one eigenmode appears to be a tremendous task. 
Selective excitation techniques as holographic masks of the modes or intracavity filtering are 
not appropriate in our case due to the complex patterns of the modes of the chaotic fiber. An 
other technique of selective excitation consists in controlling the transverse wavenumber that 
propagates along the fiber. This can be done by tuning the angle of a plane wave illumination 
at the input of the fiber. Nevertheless, due to the finite aperture of the fiber, diffraction occurs 
leading to a broadening of the value of the transverse wavenumber that propagates along the 
fiber. For a plane wave exp(— in t x) launched along the direction of the 2-bounce PO with a 
value of n t — 23.45/i?, we evaluate the spectral broadening to = 14/ R. As a consequence, 
almost 50 modes are excited preventing any possibility to maintain the modal structure of the 
scar along the propagation. To identify the guided modes along the propagation, we calculate 
what we call a pseudo-spectrum of our system. This pseudo-spectrum is derived from a stan- 
dard semiclassical procedure developped by Heller [33] and adapted to our numerical algorithm 
[34] . The pseudo-spectrum is extracted from the Fourier transform of the correlation function 
C(z) obtained from the overlap of the initial field ipo( r ) with the propagating field ^(r, z). The 
pseudo-spectrum of a plane wave illumination exp(— iK t x) with n t = 23.45/i? together with 
the NF and FF intensity at the output of the fiber are plotted in figure ([9]). In the NF figure, 
the intensity is statistically uniformly distributed over the transverse cross section of the fiber 
showing no localisation along the 2-bounce PO associated to the expected scar mode. The 
loss of information about the direction of the initial plane wave can clearly be seen in the FF 
intensity figure where the light tends to fill a ring whose radius is centered on the input trans- 
verse wavenumber n t . These pictures are actually associated to a superposition of the roughly 
50 modes initially excited as can be seen in the pseudo spectum centered on the value of the 
transverse wavenumber associated to the scar. 

Introducing doping elements in an optical fiber is a standard process to obtain optical am- 
plification for instance. As the scar modes are characterized by an enhancement of intensity 
along some particular direction in the transverse cross section of the fiber, we introduce active 
ions along that direction to perform a selective amplification through the strong interaction of 
the field with the doped area. With the experiment in mind, we have numerically investigated 
realistic systems. What can be done experimentally is to introduce active rare earth ions in a 
circular area centered on the axis of the 2-bounce PO. To make the doping area efficient for our 

1 The actual fiber has been fabricated in our institute by W. Blanc, in collaboration with P. Oupicky 
from Plasma Institute of Physics in Turnov (Czech Rep.) 
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(a) 



t 



(b) 




Kt (1/R) 



Fig. 9. NF intensity (a) and FF intensity (b) at the output of the passive fiber for a plane wave 
illumination in the direction of the 2-bounce scar with Kt — 23.45(l/i?) and the pseudo-spectrum (c). 



purpose, two conditions have to be fulfilled: a)the location of this area has to be optimized to 
favor the amplification of the 2-bounce PO scar, b) this area should have the smallest feasible 
refractive index difference with the core refractive index to prevent light from being guided 
in the doped area. As a consequence, the doped area is located close to the position of the 
self-focal point of the 2-bounce PO that also corresponds to a maximum of intensity of the scar 
modes (see figure HU|). Experimentally, the introduction of lOOOppm of Ytterbium ions leads 
to an increase of the refractive index equal to 2.10 -4 . In spite of this small refractive index 
difference between the doped area and the D-shaped core, light is actually guided in the doped 
area and amplified. To overcome this problem, we use a Gaussian refractive index profile which 
seems to support no guided modes. 




Fig. 10. Schematic view of the transverse crosssection of the active fiber with an Ytterbium doped 
area. 
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In this kind of system, optical amplification occurs due to the transfer of energy from an optical 
pump wave at the wavelength, X p = 980nm in our case, to the signal to be amplified at the 
wavelength A s = 1020nm. Ytterbium is a four-level system that permits to use different wave- 
lengths for the pump and the signal preventing direct signal reabsorption. As a consequence, 
Ytterbium in a silica matrix can be treated as a three energy-level system with a metastable 
third energy level. 

Using the evolution equations of the population density, one can easily obtain the following 
expressions for the evolution of the pump and signal intensities [35]: 

-CTpaN^Ipiz) (7) 

(y S a(r ]s N 2 -N 1 )I s (z) (8) 

where I p and I s are the intensity of the pump and of the signal respectively, <j sa the absorption 
cross section of the signal, a pa the absorption cross section for the pump and r] = cr se /a sa with 
<j se the emission cross section of the signal. As soon as the intensity of the signal I s is low 
compared to the intensity of the pump I pi Eq. 0only depends on I p and the two equations are 
decoupled so that they can be treated separately. This case is fulfilled in our simulation where 
the initial signal power is P So = 100/iW and the pump power equals to P Po = 3W. 
The Bpm algorithm is modified to take equations 17181 into account to simulate the coupled evo- 
lutions of the pump and the signal. To make an efficient transfer of energy from the pump to 
the signal, one has to maximize the number of excited pump modes: a focused laser beam at the 
pump wavelength is therefore launched in the multimode core. A great number of modes of the 
D-shaped core are thus excited and the overlap of their fields with the doped area is responsible 
for the energy transfer from the pump to the signal. In previous works, we have shown that 
the chaotic D-shaped geometry is an efficient shape for this purpose [36J. The evolution of the 
pump along the fiber follows an early exponential decay due to absorption in the doped area 
that tends to saturate as soon as the pump power is no longer able to fulfill the population 
inversion condition rj s N2 — Ni > 0. 

We are now interested in the evolution of the signal associated to a plane wave illumination in 
the direction of the 2-bounce PO of the selected scar with n t = 23 A5/R. In figure El we plot 
the NF and FF intensities computed for a distance z max « 2m associated to the maximum gain 
value. 
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(a) (b) 

Fig. 11. NF intensity (a) and FF intensity (b) of the signal at the output of the active fiber for 
launching initial condition in the direction of the scar with Kt — 23.45/ R. 



The intensity figures clearly show that the 2-bounce PO direction is privileged: in the NF inten- 
sity, high intensities concentrate on the horizontal direction associated to the scar, and in the 
FF intensity, two symmetric spots of intensity arise, located on the value of the initial hit. The 
assumption, that the expected scar mode seems to emerge thanks to amplification, is confirmed 
by the pseudo-spectrum. The figure [12] represents the evolution of the pseudo-spectrum along 
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the propagation. It is deduced from pseudo-instantaneous spectra computed for three sample of 
propagation of length 5cm along the 4-m fiber. From a broad spectrum for z = Om, we observe 
that a selection of few peaks tends to emerge. Only some modes are being amplified along 
the propagation, and the mode associated to the scar with value of K t = 23 A5/R distinctly 
dominates the spectrum. 
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(a) (b) 
Fig. 12. Evolution of the intensity spectrum along the propagation (a) at z=0cm, (b) z=80cm 
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Fig. 13. Spectrum calculated around the maximal value of the gain and the calculated modes associated 
to its main peaks. 



To clearly identify the excited modes, we perform a filtering on the FF by applying a ring with 
mean radius n t and a Gaussian profile, K t being associated to one of the main peaks of the 
spectrum. The corresponding NF is deduced from the inverse Fourier transform. By applying 
the previous procedure to the value of n t related to the scar of order p=10, we have checked 
that the spatial distribution of the intensity is associated to the corresponding scar mode (fig 
IT3]) . We have also checked that the peak labelled by n t = 10.9/ii corresponds to the intensity 
pattern of the scar of order p=4 along the 2-bounce PO. The NF intensities obtained by apply- 
ing the filtering procedure on different values of n t respectively associated to the scar modes 
of order p=A and p=10 are plotted in an insert on figure [13] and compared to the correspond- 
ing modes calculated for a waveguide with Dirichlet conditions. The intensity patterns for the 
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two filtered-NF are similar to the individual calculated scars: the scar modes are amplified to 
detriment of the ergodic modes. 

To check the critical effect of the location of the gain region on scar amplification, we have 
changed the location of the doped area. By moving the location of the doped region away from 
the self- focal point of the 2-bounce PO, we observe no scar amplification confirming that the 
gain region must have a maximal overlap with the scar to be efficient. 

The localization of light along the 2-bounce PO clearly demonstrates that the scar modes 
have been excited and enhanced thanks to the interaction with the doped region. 



4 Conclusion 

The optical fiber with a truncated cross section appears to be a good candidate to image experi- 
mentally wave chaos: complementary informations about the field that propagates in the chaotic 
billiard are deduced from the NF and FF intensity. The ergodic behavior of generic modes has 
been proved experimentally and scar modes have first been observed in an optical fiber. Due 
to the inescapable diffraction at the fiber input, the properties of individual scar modes cannot 
be studied experimentally. Nevertheless, their signature on the probability distribution of the 
intensity appears to be pregnant: our preliminary investigations have shown that deviations to 
the ergodic theory are to be expected as soon as a scar is excited on a superposition of modes. 
To perform a selective excitation of scar modes in our multimode structure, we have investi- 
gated the influence of a localized gain area on the propagation of a plane-wave illumination. 
Our numerical results clearly show that scar modes are selectively amplified. The experimental 
validation is in progress and will be the subject of a next publication. 
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